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Abstract 

We present a method for calculating the full current noise spectrum S(u>) for the class of nano-electromechanical systems 
(NEMS) that can be described by a Markovian generalized master equation. As a specific example we apply the method to 
a quantum shuttle. The noise spectrum of the shuttle has peaks at integer multiples of the mechanical frequency, which is 
slightly renormalized. The renormalization explains a previously observed small deviation of the shuttle current compared to 
the expected value given by the product of the natural mechanical frequency and the electron charge. For a certain parameter 
range the quantum shuttle exhibits a coexistence regime, where the charges are transported by two different mechanisms: 
shuttling and sequential tunneling. In our previous studies we showed that characteristic features in the zero-frequency noise 
could be quantitatively understood as a slow switching process between the two current channels, and the present study 
shows that this interpretation holds also qualitatively at finite frequency. 
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1. Introduction 

A decade of advances in microfabrication technol- 
ogy has pushed the typical length scales of electrome- 
chanical systems to the limit, where quantum mechan- 
ical effects of the mechanical motion must be taken 
into account [1] . Such nano-electromechanical systems 
(NEMS) exhibit a strong interplay between mechan- 
ical and electronic (or magnetic) degrees of freedom, 
and their electronic transport properties reflect this in- 
terplay in an intricate manner. 
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A modern trend in transport studies of mesoscopic 
systems has been to not only consider the current- 
voltage characteristics of a given device, but also to 
examine the noise properties, or even the higher cu- 
mulants (i.e. the full counting statistics (FCS)) of the 
current distribution [2,3]. The current noise, either its 
zero- frequency component or the whole frequency spec- 
trum, provides more information than just the mean 
current and can be used to discern among different pos- 
sible mechanisms resulting in the same mean current. 
While the noise spectra in generic mesoscopic systems 
have been studied well over a decade, it is only very 
recently that the study of noise spectra of NEMS has 
been initiated [4,5,6]. 
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The aforementioned three studies deal with the noise 
spectra of a classical shuttle, a classical nanomechan- 
ical resonator coupled to a single electron transistor 
(SET), and the Ceo-SET in a strong electromechani- 
cal coupling regime, respectively. The first two stud- 
ies [4,5] found peaks in the current noise spectra at 
the first two multiples of the mechanical frequency. For 
low bias voltages and strong electromechnical coupling 
the third study [6] found a power-law frequency de- 
pendence of the noise spectrum attributed to scale-free 
avalanche charge transfer processes. In all three cases 
the noise spectra revealed interesting details about the 
systems. From the technical point of view, two of the 
studies [4,6] used Monte Carlo simulations, whereas [5] 
used a model-specific numerical evaluation of the Mac- 
Donald formula (see below) . 

In this work, we present a study of the full frequency 
spectrum of the current noise of a quantum shuttle 
[7,8,9]. We extend the general formalism developed for 
the zero-frequency noise [10,11] and the FCS [12] cal- 
culations for NEMS described by a Markovian general- 
ized master equation. The presented formalism applies 
not only to the shuttle studied here but could equally 
well be used for all three systems from the previous 
studies [4,5,6] for the determination of the noise spec- 
tra. 

We apply the developed theory to compute numer- 
ically the noise spectrum of the shuttle in the deep 
quantum regime 2 as function of the damping coeffi- 
cient. The spectrum has peaks at integer multiples of 
the slightly renormalized mechanical frequency. The 
renormalization of the bare oscillator frequency as read 
off from the current spectrum explains a small but ob- 
servable deviation from the expected value of the cur- 
rent in the shuttling regime i s hut = euio /2n [8] . It turns 
out that it is the renormalized oscillator frequency u>o 
which should enter this relation. Finally, we focus on 
the low-frequency part of the spectrum when approach- 
ing the semi-classical regime for intermediate values 
of the damping, i.e. in the coexistence regime, where 
both shuttling and tunneling are effective. We use the 
frequency dependence of the spectrum for w < wo to 
identify additional qualitative evidence for the bistable 
behavior of the shuttle in this regime described by a 



2 In this regime A ~ xq, where A is the tunneling length 
(see also Eq.(l)) and xo = h/rriLUQ sets the length scale 
for the quantum mechanical zero-point motion. 




x 

Fig. 1. The quantum shuttle consists of a nanosized grain 
moving in a harmonic potential between two leads. A high 
bias between the leads drives electrons through the grain. 



simple analytical theory of a slow switching between 
two current channels (compare with Refs. [12,13]). 



2. Model 

We consider the model of a quantum shuttle de- 
scribed in [7,8,9,10]. The shuttle consists of a mechan- 
ically oscillating nanoscale grain situated between two 
leads (see Fig. 1). In the strong Coulomb blockade 
regime the grain can be treated as having a single elec- 
tronic level only. A high bias between the leads drives 
electrons through the grain and exerts an electrostatic 
force on the grain, when charged. The grain is assumed 
to move in a harmonic potential, and the oscillations 
of the grain are treated fully quantum mechanically. 
Damping of the oscillations is described by interactions 
with a surrounding heat bath. 

From the Hamiltonian of the model one can derive 
a generalized master equation (GME) resolved with 
respect to the number of electrons n that have been 
collected in the right lead during the time span to 
t. The n-resolved GME describes the time evolution 
of the n-resolved system density matrix p^ n \i), where 
the 'system' consists of the electronic level of the grain 
and the quantized oscillations. In the following we only 
need the n-resolved GME for the part of /r ' (t) that is 
diagonal in the electronic components, which reads [10] 



2 



j$ (t) =^[Hosc, f$ (t)] + (t) 

fi?(t) =^[Ho SC - eEx,p£>(t)] + C d ^ p p^(t) 

(1) 

with n = 0,1,2,... and p\[ X ^ {€) = 0. The commu- 
tators describe the coherent evolution of the charged 
(M? = <l|P (n) |l» or empty (p$ = (0|p<"'|0» shuttle 
with mass m and natural frequency lo . The electric 
field 3 between the leads is denoted E. The terms pro- 
portional to ^l(r) describe transfer processes from the 
left (to the right) lead with hopping amplitudes that 
depend exponentially on the ratio between position x 
and the electron tunneling length A. The mechanical 
damping of the oscillator is described by the damp- 
ing kernel (here at zero temperature) £ d 

ffi}] ~ ^[x,[&,P% ) ]]J = 0,i [8,10]. 

The n-resolved GME can be recast into the compact 
form [11] 

fiW = (C-l R )p (n) +l R p ( "- 1 \ P ( - 1) =0, (2) 

where we have introduced the Liouvillean £, describ- 
ing the evolution of the system density matrix p(t) = 
~Y2 n p( n ) (t), i.e. pit) = Cp(t), and the superoperator for 
the tunnel current through the right junction (taking 
e = 1), defined by its action on the density operator 

Xfl/3 = r fl e*|OXl|/3|lXO|e*. (3) 

Assuming that the system tends exponentially to a 
stationary state p stat the Liouvillean has a single eigen- 
value equal to zero with p stat being the (unique and nor- 
malized) right eigenvector which we denote by |0)) [11]. 
The corresponding left eigenvector is the identity oper- 
ator 1 which we denote by ((0|, and from the definition 
of the inner product 4 we have «0|0» = Tr(i t p stat ) = 



3 In order to obtain a Markovian GME it is necessary to 
assume that the bias applied between the leads is the high- 
est energy scale of the system [8]. Nevertheless, we consider 
the electric field between the leads as a free parameter of 
the model. 

4 We define the inner product of two supcrvcctors as 
((a\b)) = Tr(AtB) with the identification \o)) <-> 6, where 
O is a quantum mechanical operator and |o)) is the corre- 
sponding supcrvector. 



1. In terms of Xr the average tunnel current in the sta- 
tionary state can be expressed as 

(J fl )=Tr(J fl p stat ) = «6|Xi ? |0». (4) 

We define the projectors V = |0»«0| and Q = 1 - V 
obeying the relations VC = CV = and QCQ = C. In 
terms of the two projectors we can express the resolvent 
of the Liouvillean Q(— iuj) = (— iu) — £) _1 as 

G{-iu) = -—V-Q—^—;Q = -—V-11{lu), (5) 

where we have introduced the frequency dependent su- 
peroperator 1Z(uj), which is well-defined even for to = 
0, since the inversion in that case is performed only in 
the subspace where C is regular. 



3. Theory 

We consider the current autocorrelation function de- 
fined as 

C II (t',t") = ±{{Ai(t'),Ai(t")}), (6) 

where AI(t) = I(t) — (I(t)). In the stationary state 
Cn(t' , t") can only be a function of the time difference 
t = t' — t" , and we thus write 

Cu(t) = l({Ai(t),Ai(0)}). (7) 

The current noise spectrum is the Fourier transform of 
Cu(t), i.e. 

/OC 
dtCu{t)e iuJt . (8) 
- OC 

In order to calculate the current noise measurable in, 
say, the right lead one must recognize that the cur- 
rent running in the lead is a sum of two contributions, 
namely the tunnel current through the right junction 
and a displacement current induced by electrons tun- 
neling between leads and grain. This is reflected in the 
Ramo-Shockley theorem [2] 

i = c L i R + c R i L . (9) 

Here / is the current operator for the current running 
in the lead, whereas Il(r) is the operator for the tunnel 
current through the left (right) junction, and cum is 
the relative capacitance of the left (right) junction in 
the sense cl + cr = 1. Combining the Ramo-Shockley 
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theorem with charge conservation leads to an expres- 
sion for the current noise measured in the lead [14] 

Sii(uj) = ClSi r i r (lu) + crSi l i l (lu) — c l crL) 2 Snn{u), 

(10) 

where N = |1)(1| is the occupation number operator 
of the electronic level of the grain. In the following 
we neglect any dependence of the capacitances on the 
position of the grain and consider the symmetric case 

Cl = Cr = i. 

The two first terms of Eq. (10) can be evaluated 
using the methods developed by MacDonald[15]. The 
starting point of the derivation is the property Cn(t) = 
Cn(—t), which immediately leads to 

{■ oo 

S lRlR {u)= dtd^me^ + e-^). (11) 
Jo 

Let us consider the first term 



and consequently 

Si r i r (cu) = / dt — (AQ R (t))(cu sin cut+e cos ut) 
Jo dt 

We now make use of the fact that 



e 

(18) 



(AQ 2 R (t)) = {Ql(t)) - (Q R (t)) 2 = (n 2 {t)) - (n(t)) 2 , 

(19) 

with (n a (t)) = J2^ n a P n (t),a = 1,2, where P n (t) 
is the probability of having collected n = 0, 1, 2, . . . 
electrons in the right lead during the time span to t. 
This finally leads us to the commonly used form of the 
MacDonald formula (cf. [5,14,15]) 

S lRlR {u>) = w j dtsin(uot)- [<n 2 (t)> - (n(t)} 2 ] , 

(20) 

where the regularization 



'IrIr 



: / dtCi R i R (t)e iw \ (12) lj sin(wt) -> ± ((w + ie)e i{u+ie)t - [u - ie)^ 1 ^'*) 

Jo 2i 



the second term, Sj j (w), follows analogously. Defin- 
ing Q R (t) as the operator of charge collected in the 
right lead in the time span to t we have 

AQ R (t) = Q R (t) - (Q R (t)) = / dt'AI R {t'), (13) 



= (lj sin cut + e cos uit)e s ,e^0 



(21) 



and we can express the current autocorrelation func- 
tion as 

Introducing the convergence factor s —> + and per- 
forming the integration by parts in Eq. (12) we get 



is implied. Only the proper treatment of the regular- 
ization ensures correct results including the to = case 
where the zero-frequency noise MacDonald formula is 
recovered (by using the Laplace transform identity for 
e^0+) 



[■OO J 

Si R i R (0) = e I dte~ st f t [(n 2 (t)) - (n(t)f 



d 



(22) 



= Tt [<" (*)) " ("(*)> 



s tin^)= / dt{{AQ R (t),AI R (0)}) 
Jo 



x ^ + i(u+ie)t 

2i 



(15) 



dt " J lt-»oo 
In order to evaluate the current noise spectrum we 
now introduce the quantity 



Jo 



dte 1 ' 



±{n 2 (t))-2{n(t))±{n(t)) 



Since {{AQ R (t),AI R (0)}} = ({AQ R (t), AI R (t)}) = 



-^(AQ R (t)) in the stationary state 5 we can write 



Similarly, we find 

sj Rl M = -f 

Jo 



(23) 

with either S(uj) = ImS(w) or S{u) = (S(lu) + 
S(—lo)) /2i and evaluate it along the lines of [11]. Since 
P„(t) = Tr(p (n) (t)), Eq. (2) leads to (keeping in mind 
that Tr(£») = 0) 

P n (t)=Tr[l R (p^- 1 \t)- 



(*))], 



dt 



2i 



(17) 



5 The first equality follows from a simple substitu- 
tion t — > t — t in the chain ({AQ R (t), AI R (0)}) = 
J*dr({Ai R (r),Ai R (0)}) = f* dr({AI R (t), AI R (r)}) = 
{{AQ R (t),Ai R (t)}). 



and as shown in [11] 

±(n(t)) = T±(X R p(t)) = {(0P*\0)), 



(24) 



(25) 



f t {n 2 (t)) = 2Tr{lRY,np (n \t)] + fe|0», (26) 
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where we have used p(t) = p stat , since we are consider- 
ing the stationary limit. The sum entering Eq. (26) is 
evaluated by introducing an operator-valued generat- 
ing function defined as F(t, z) = £^=0 P (n) (t)z n , from 
which we get 

d 



dz 



F(t,z)\ =J2np (n) (t). 

I 2—1 A * 



(27) 



From the definition of the Laplace transform, F(s, z) = 
/ °° dtF(t, z)e~ st , we see that the integration in Eq. 
(23) can be considered as a Laplace transform evalu- 
ated at s = — iui + e (remember the proper regular- 
ization; from now on we skip explicitly mentioning the 
e-factors). In [11] it was shown that 

£-F(s=-iu,z)\ 



=g(-iw)i R g(-iu)p(o) 

+ g(-ico)J2np (n) (0). 



(28) 



Again, we have p(0) = p stat , and moreover we as- 
sume the factorized initial condition [14,16] p^ n '(0) = 
<5onp stat , i-e. we start counting the charge passing 
through the right junction at t — and the system is 
in its stationary state. Now, combining Eqs. (5, 23-28) 
and having in mind that -Pp stat = p stat , Qp stat = 0, 
straightforward algebra leads to 



S(u) = i («0|X A |0» - 2({b\T R TZ(uj)l R \0)}) . (29) 
We thus arrive at 

SwnH = «o|Xb|o» - 2R C [{{dp R n{ u )x R \o))] 



«0|Jb|0» - 2«0|X fl 



c 



C 2 + u> 2 



(30) 
(31) 



For the left junction one similarly finds 

Si l i l {u) = «6|Xl|0» - 2Rc [{{0\1 l 11(lu)1 l \0): 
with 

Ji/3 = rie-*|l)(0|p|0Xl|e-*. (32) 

For the evaluation of the charge-charge correlation 
function Snn(u) we note that N is a system operator, 
and the quantum regression theorem thus applies [17]. 
Following [11] we immediately get 

Snn{co) = -2R C [{(b\Afn(w)Ar\o))] , (33) 

having introduced the superoperator Af corresponding 
to N, defined as 



Collecting all terms in Eq. (10) we finally obtain the 
expression for the current noise measured in the leads 
for the symmetric setup (cl = c R — i) 

S n (uj) =«6|Xk|0» + ^Rc [«0^(w)7V|0»] 

- Rc [{{0\1 r 1Z(u;)1 r +1 l TZ(uj)1l\0})] . 

(35) 

We notice that for uj = we get the previous result 
[10,11,12] 

5/7 (0) = «0|X^|0» - 2((0|2: fl ^(0)Jfl|0», (36) 

since the zero-frequency tunnel current noise is the 
same at both junctions and the second term is real [11]. 

The numerical evaluation of Eq. (35) is only possi- 
ble by truncating the number of oscillator states. As 
in previous studies [8,10,12] we retain the 100 lowest 
oscillator states, which however still leaves us with the 
task of dealing numerically with the matrix represen- 
tations of the relevant superoperators, which are of size 
20000 x 20000. As explained in [11] the stationary den- 
sity matrix p stat (or j0))) can be found using the Arnoldi 
iteration scheme, and 1Z(lo) can be evaluated using the 
generalized minimum residual method (GMRes). Both 
methods are iterative and rely crucially on an appro- 
priate choice of preconditioner to ensure that the it- 
erations converge and to speed up the computation. 
It should be stressed that finding a suitable precondi- 
tioner for a given problem is by no means simple. For 
finding p stat and 11(0) it turns out that the Sylvester 
part of the Liouvillean, which is the part that can be 
written C sy \ v p = Ap + pA* , is well-suited for precon- 
ditioning [11]. This preconditioner separates the zero 
eigenvalue from the rest of the spectrum of C leading 
to a decrease in computation time. 

For finding 1Z(lu) the original preconditioner must be 
modified in order to separate the relevant eigenvalue 
from the rest of the spectrum. A reasonable choice is 
the superoperator M defined as 



M P ={A+ —)p + p (A-—y. 



(37) 



Jvp=|l)<l|p|l)<l|. 



(34) 



For the range of parameters discussed in the present 
paper this choice of preconditioner was sufficient for 
convergence, however the obtained computational 
speedup is considerably smaller than the speedup 
provided by £ sy i v , when computing 1Z(0). This obser- 
vation combined with the fact that GMRes fails to 
converge for certain parameters in the semi-classical 
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Fig. 2. The ratio between the current noise and the current 
F(w) = Su(w) / (I) as function of the damping 7 and fre- 
quency u. The other parameters arc = T^j = 0.05luq, 
A = xq, d = eE/mu>Q = 0.5xq, where xq = y/ fi/muo . 
Peaks are seen at ui ~ 1.03o;o, 2.06u;o, 3.09<^o- The peak at 
u) ~ 2.06^0 reaches values of F(w) ~ 20 (not shown) for 
7 = 0.02 and decreases monotonously with increasing 7 to 
-F(a;) ~ 6 for 7 = 0.09. The insert shows a representative 
curve (7 = 0.05). 

regime indicates that the identification of the optimal 
preconditioner for the problem at hand remains an 
open problem. 



4. Results 

In Fig. 2 we show results for the current noise spec- 
trum of the quantum shuttle in the deep quantum 
regime. In accordance with previous studies [4,5] we 
find peaks at integer multiples of the mechanical fre- 
quency. A close look at the spectrum reveals a slight 
renormalization of the mechanical frequency with the 
peaks appearing at ui — (Do , 2(Do , 3d>o , where ujq ~ 
1.03ojo- In the shuttling regime the current is expected 
to saturate at a value expressed as one electron per 
cycle of the mechanical vibrations. For a shuttle with 
mechanical frequency u>o this implies that the satu- 
rated shuttle current is /shut = UJo/lm ~ 0.159cjo- For 
the given parameters the numerical calculation yields 
a slightly higher value, namely /shut = 0.164cjo = 
0.159 x 1.03wo — luo/2h, and this can now be under- 
stood in light of the observed renormalization of the 
mechanical frequency. 

In [12,13] it was shown that an observed giant en- 
hancement of the zero-frequency noise [10] in the co- 



existence regime of a shuttle approaching the semi- 
classical regime can be understood in terms of a simple 
model of a bistable system switching slowly between 
two current channels (shuttling and tunneling) . Denot- 
ing the currents of the two channels 6 as 7s and It, 
respectively, and the switching rates Fs<_t and Ft<— s, 
one can show (following [13]) that the ratio between 
the current noise spectrum and the current (in the 
zero- frequency limit known as the Fano factor) F(ui) = 
S(lj) 1 1 for the bistable system has the Lorentzian form 



F(co) = - 



2 r £ 



(Is - It) 2 



where 



/ Fs^t + F T —s (Fs<-t + r T *-s) 2 + 1 



IsFs^t + Ir^T^i 



(38) 
(39) 



Fs^t + Ft^s 
The two switching rates, Fs^t and Ft^s, can be ex- 
tracted from the numerical values of current and zero- 
frequency noise, and by comparing Eq. (38) with the 
noise spectrum obtained numerically, one can perform 
another independent test of the hypothesis that the 
shuttle behaves as a bistable system in the coexistence 
regime. 

In Fig. 3 we show numerical results for the low- 
frequency current noise of the quantum shuttle in the 
coexistence regime. Together with the numerical re- 
sults we show the analytic expression for the current 
noise of the bistable system (Eq. (38)) with rates ex- 
tracted from the numerical values of current and zero- 
frequency noise. It should be noted that the agreement 
between the numerical and analytic results could only 
be obtained by assuming that the shuttling current for 
the given values of the damping is not fully saturated 
to the value /shut = £>o/27r. The current noise spectrum 
thus provides us with qualitative evidence for the shut- 
tle behaving as a bistable system in the coexistence 
regime, while it, however, leaves us with an open ques- 
tion concerning the saturation of the shuttling current. 



5. Conclusion 

We have presented a theory for the calculation of 
the current noise spectrum of a large class of nano- 



, and It 



r r .r f 



with Fr 



2tt> " r L +r R ' 
r R e h / muJ 0* 2 e 2eB/mwgA j j^. = r ie V™"(l* 2 [10,13]. 
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Fig. 3. The ratio between the current noise and 
the current F(u>) = Sn(io)/(I) for low frequen- 
cies (a; <C a;o). The parameters are 7 = 0.035u;o 
(lowest curve), 0.04ojo, 0.045u;o, 0.05wo (topmost curve), 
T L = F R = O.OIijjo, A = 1.5xo, d = eE/mu}Q = 0.5a;o, 
where xq = h/mu>o . The circles indicate numerical re- 
sults, while the full lines indicate the analytic results for 
the current noise spectrum of a bistable system. It should 
be noted that in order to obtain the agreement between 
the numerical and analytic results it is necessary to as- 
sume that the shuttling current for the given values of 
7 is not fully saturated to the value J s h u t = 1.03a>o/27r. 
Corresponding to the different values of 7 we have used 
Ishut = 1.01oj /27r(for 7 = 0.03cj ), 1.00wo/2tt, 0.98w /2tt, 
0.94uo/27r(for 7 = 0.05c<Jo); respectively. 



electromechanical systems, namely those that can be 
described by a Markovian generalized master equa- 
tion. As a specific example we have applied the the- 
ory to a quantum shuttle. For the quantum shuttle 
numerical calculations of the current noise spectrum 
in the deep quantum regime revealed a slight renor- 
malization of the mechanical frequency - this, in turn, 
explains an observed small deviation of the shuttle 
current compared to the expected value given by the 
product of the natural mechanical frequency and the 
electron charge. When approaching the semi-classical 
regime the low-frequency noise served as another evi- 
dence for the quantum shuttle behaving as a bistable 
system, switching slowly between two current channels, 
thus supporting this claim, previously based on the cal- 
culation of the full counting statistics of the quantum 
shuttle. The theory presented here has a broad range 
of applicability, encompassing the few previous studies 
of the current noise spectra of nano-electromechanical 
systems. 
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